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ABSTRACT 

We present photometry of four transits of the planetary system HAT-P-13, obtained us- 
ing defocussed telescopes. We analyse these, plus nine datasets from the literature, in or- 
der to determine the physical properties of the system. The mass and radius of the star 
are Ma = 1.320 ± 0.048 ± 0.039 Mq and Ra = 1.756 ± 0.043 ± 0.017 R© (statistical 
and systematic errorbars). We find the equivalent quantities for the transiting planet to be 
Mb = 0.906 ± 0.024 ± 0.018 Mjup and i?b = 1.487 ± 0.038 ± 0.015 Rjup, with an equi- 
librium temperature of T^^ — 1725 ± 31 K. Compared to previous results, which were based 
on much sparser photometric data, we find the star to be more massive and evolved, and the 
planet to be larger, hotter and more rarefied. The properties of the planet are not matched by 
standard models of irradiated gas giants. Its large radius anomaly is in line with the observa- 
tion that the hottest planets are the most inflated, but at odds with the suggestion of inverse 
proportionality to the [^] of the parent star. We assemble all available times of transit mid- 
point and determine a new linear ephemeris. Previous findings of transit timing variations in 
the HAT-P-13 system are shown to disagree with these measurements, and can be attributed 
to small-number statistics. 
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1 INTRODUCTION 

The discovery of the HAT-P- 1 3 bv' Sakos et all 1200?) elicited sub- 
stantial interest. The transiting extrasolar planet (TEP) HAT-P-13 b 
and its host star HAT-P-13 A are rather typical examples of these 
objects, with two exceptions. Firstly, the star i s the second-mos t 
metal-rich known to host a TEP after XO-2 l lBurke et al.ll2007l) 
(but see the discussion on host star [|^] values in Enoch et al. 
I2OIII) . Secondly, there is a third component in the system which 
is clearly detected in the radi al velocity measurements (RVs) of the 
host star teakos et alj |2009h . HAT-P-13 c has an orbit with a pe- 
riod of 446.22 ± 0.27 d, an eccentricity of 0.6616 ± 0.005 2 and 
a minimum mass of 14.28 ± 0.28 Mjup Iwinn et alj|2010bl) . This 
object is expected to induce transit timing variations (TTVs) within 
the HAT-P-13 A,b syst em, which potentially allow the structure o f 
the planet to be probed i Mardling & Linl2004lBatvgin et al.l2"009h . 
HAT-P-13 is unfortunately not the best system for such analyses, as 
its relatively long and shallow transits are not conducive to precise 
timing measurements. 

An observati on of the Rossiter-McLaughlin effect by 
IWinn et al.l j2010bh has shown that the projected angle between 
the orbital axis of HAT-P-13 b and the rotational axis of the par- 



See: |http : / /www, astro . keele . ac.uk/~ jkt/tepcat/| 



ent star is consistent with zero. This is in line with the fact that 
misaligned axes are only fo und for TEP systems containing a star 



misaligned axes are only to und tor 1 hi- systems containing a star 
hotter than roughly 6 250 K jWinn et id]|2010al : ISchlaufmaijboid : 
lAlbrecht et"al]|201 ih . IWinn et al.l (120101^ also found a long-term 
drift in the radial velocity measurements of the star, which may be 
the signature of a fourth component to the system on an orbit of 
much longer period. 



ISzabo et alj ( I2OI0I) attempted to detect a transit of the third 
body, at a time predicted by I Winn et al. I ( l2010bh , but were not 
successful. Their observational campaign included scrutiny of two 
transits of HAT-P- 13 b, whose times of occurrence agreed well 
with the predicted timings. IPaletal .1 ( 1201 l|) subsequently presented 
photometry of three transits, all of whose midpoints fell earlier 
than ex pected according t o an ephemeris b uilt on the observa- 
tions of iBakos et al] ( |2009|) and lSzabo etakl Uo iO). Pal et al. in- 
terpreted this as evidence of TTVs. lNascimbeniet aL C201 Ih) pre- 
sented high-speed photometry of five transits in early 2011 , taken 
as part of the TASTE project ^Nascim beni et alj I2OI lah . They 
confirmed that a linear ephemeris could not explain the avail- 
able transit timings, and postulated that a sinusoidal TTV with 
an amplitude of 0.005 d was a good match to the available transit 
ti mings. Sinusoidal TTVs have previously been seen for WASP- 
3 /Pollacco et al. 2008; Maciejewski et al. 2010) and WASP-10 
( ichristian et al., ,2009. : .Maciejewski et aL .2011i) but have not yet 
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Table 1. Log of the observations presented in tliis work. A^obs is the number of observations and 'Moon ilium.' is the fractional illumination of the Moon at 
the midpoint of the transit. The aperture sizes are the radii of the software apertures for the object, inner sky and outer sky, respectively. 



Transit 


Date 


Start time 


End time 




Exposure 


Filter 


Airmass 


Moon 


Aperture 


Scatter 






(UT) 


(UT) 




time (s) 






ilium. 


sizes (px) 


(mmag) 


Cassini 


2011 02 06 


19:21 


00:37 


127 


120 


Thuan-Gunn i 


1.17^ 1.00^ 1.10 


0.129 


25, 35, 55 


0.77 


Cassini 


2011 04 17 


18:58 


23:50 


147 


90 


Thuan-Gunn i 


1.03 2.11 


0.998 


18, 30, 50 


0.93 


Portalegre 


2011 01 31 


22:13 


05:21 


128 


180 


Cousins / 


1.11 -s> 1.01 1.64 


0.041 


10, 35, 45 


2.47 


Portalegre 


2011 02 03 


20:53 


02:55 


119 


150 


Cousins / 


1.25 1.01 1.16 


0.008 


10, 35, 45 


2.53 



Table 2. Excerpts of the light curves of HAT-P-13. The full dataset will be 
made available at the CDS. 



.00 



Telescope 


BJD(TDB) 


Diff. mag. 


Uncertainty 


Cassini 


2455599.30664 


0.00061 


0.00090 


Cassini 


2455599.52535 


-0.00085 


0.00078 


Cassini 


2455669.29046 


-0.00014 


0.00108 


Cassini 


2455669.49355 


0.00066 


0.00117 


Portalegre 


2455596.37076 


-0.5552 


0.0018 


Portalegre 


2455596.62212 


-0.5460 


0.0017 


Portalegre 


2455593.42575 


-2.0886 


0.0022 


Portalegre 


2455593.72313 


-2.0970 


0.0024 



been confirm ejfj The putative T TVs for HAT-P-13 have been 
challenged bv iFulton et al.l ( l201lh . who presented observations of 
ten transits over two observing seasons. They found that a lin- 
ear ephemeris was an acceptable match to all transit timing mea- 
surements, with th e exception of the first of the two obtained by 
ISzab6etal.l ( l2O10h . 

Th e physical propertie s of t he HAT-P- 13 system have been de- 
rived bv lBakos et^ J2009l) and lWinn et al] i2010bl) . who used the 
same light curves. Since these studies a wealth of new photometric 
data has been gathered. This has been used to investigate putative 
TTVs, but has not been brought to bear on improving the physical 
properties of the system. In this work we present new photometry 
covering four transits and use all available high-quality photometry 
to measure refined physical properties of HAT-P-13. 



2 OBSERVATIONS AND DATA REDUCTION 

Two full transits of HAT-P-13 were observed with the BFOSC im- 
ager mounted on the 1.52 m G. D. Cassini Telescop^fl at Loiano 
Observatory, Italy. We used a Gunn i filter and autoguided through- 
out. We had to reject a small number of datapoints in both transits, 
as they were affected by pointing jumps which compromised the 
data quality. A summary of our observations is given in Table[T]and 
the full data can be found in Table|2] 

The telescope was defocussed so the point spread functions 
(PSFs) resembled annuli of widths 15-25 pixels, in order to reduce 
the light from the target and comparison stars to a maximum of 
roughly 35 000 counts per pixel. This approach reduces the sus- 
ceptibility of the data to flat-fielding noise and increases the ef- 
ficiency of the observations. A detailed description of the defo- 



2 A sinusoidal TTV was claimed for OGLE-TR-1 1 1 bv lDi'az et all j2008l) 

but has been refuted by Adams et al. (2010) 

^ Information on the 1.52m Cassini Telescope and BFOSC can be found 
at jhttp : / / www .bo.astro.it/loiano/. 
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Figure 1. New data presented in this work, compared to the best JKTEBOP 
fits using the quadratic LD law. The dates of the Ught curves are labelled 
using the format month/day. The residuals of the fits are plotted in the lower 
half of the figure, offset from zero. 



cussing method can be found in ISouthworth eTaP (l2009allbl), and 
an inst ance of its use with the Cassini telescope in Southworth et al.l 
( l2010h . Several images were taken with the telescope properly fo- 
cussed, and used to verify that there were no faint stars within the 
defocussed PSF of HAT-P-13. 

Data reduction was undertaken using standard methods per- 
taining to aperture photometry. Software aperture positions were 
specified by hand but shifted to account for pointing variations, 
which were found by cross-correlating each image against the ref- 
erence image used to place the apertures. We found that the re- 
sults were insensitive to the choice of aperture sizes (within reason) 
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Figure 2. Phased light curves of HAT-P-13 compared to the best JKTEBOP fits using the quadratic LD law (left panel). They are shown in the same order as in 
Table[3] The residuals of the fits are plotted in the right panel, offset to bring them into the same relative position as the corresponding best fit in the left panel. 



and to whether flat fields were used in the data reduction process. 
Differential-photometry light curves were obtained with respect to 
an opti mal ensemble of four com parison stars constructed as out- 
lined bv lSouthworth et ^ i2009ah . The times of observation were 
converted to barycent ric Julian date on the TDB timescale, using 
the IDL procedures of lEastman et alj j2010l) . 

Two full transits were observed by JG from CROW- 
Portalegre, Portugal, using an f/5.6 30 cm Schmidt-Cassegrain tele- 
scope, a KAF1603 CCD camera operating at a plate scale of 
1.12" px^^, and a Cousins / filter. The data were reduced by stan- 
dard methods, using median-combined bias, dark and flat-field cal- 



ibration observations. Aperture photometry was performed with C- 
Munipacl0 and differential-magnitude light curves obtained with 
respect to an ensemble of five comparison stars. 



3 LIGHT CURVE ANALYSIS 

We have analysed the Cassini and literature photometric observa- 
tions of HAT-P-13 by the methods of the Homogeneous Studies 

* |http : //c-munipack . sourcef orge ■ net/| 
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Table 3. Parameters of the JKTEBOP fits to the light curves of HAT-P-13. The final parameters correspond to the weighted mean of the results for the ten light 
curves. 



Sourcs 




k 


i ("I 






Cassini -i-band 


0.211 ±0.014 


0.0932 ±0.0015 


81.6± 1.0 


0.193 ±0.013 


0.0180 ±0.0014 


Bakos rLWO ?-band 


0.193 ±0.01 1 


0.08592 ± 0.00084 


82.61 ± 0.87 


0.178 ± 0.010 


0.0153 ± 0.0010 


Szabo V-band 


0.192 ±0.030 


0.0836 ± 0.0046 


83.8 ±3.2 


0.178 ±0.028 


0.0148 ±0.0031 


Szabo ij-band 


0.207 ± 0.024 


0.0872 ± 0.0026 


81.7± 1.8 


0.190 ±0.022 


0.0165 ±0.0019 


Pal _R-band 


0.217 ±0.034 


0.0895 ± 0.0073 


81.5±2.6 


0.200 ±0.031 


0.0178 ±0.0036 


Pal /-band 


0.210±0.019 


0.0887 ± 0.0032 


81.5± 1.3 


0.193 ±0.017 


0.0171 ±0.0018 


Nascimbeni _R-band 


0.2103 ±0.0048 


0.08668 ± 0.00052 


81.97 ±0.36 


0.1852 ±0.0044 


0.01606 ±0.00042 


Fulton FTN Z-band 


0.231 ±0.023 


0.0882 ±0.0041 


80.2 ± 1.6 


0.212 ±0.021 


0.0187 ±0.0021 


Fulton FLWO «-band 


0.208 ±0.016 


0.0868 ± 0.0023 


81.5± 1.2 


0.191 ±0.015 


0.0166 ±0.0017 


Fulton Sedgwick i-band 


0.199±0.017 


0.0873 ± 0.0029 


82.4 ± 1.4 


0.183 ±0.016 


0.0159 ±0.0017 


Final results 






81.93 ± 0.26 


0.1863 ± 0.0034 


0.01622 ±0.00034 


_Bakos et al. (2009) 


0.1856 


0.0844 ±0.0013 


83.4 ±0.6 


0.1712 ±0.0076 


0.01445 


Winn et al. (2010b) 


0.1839 


0.08389 ±0.00081 


83.40 ±0.68 


0.1697 ±0.0072 


0.01424 


Fulton etal. (2011) 


0.1967 


0.0855 ±0.0011 


82.45 ± 0.46 


0.1812 ±0.0056 


0.01549 



project ( ISouthwortlJ[200l I2009L I2OI0I . l201lh . which are briefly 
summarised below. The light curves and their best-fitting models 
are shown in Fig.[T]for the data presented in this paper, and in Fig.|2] 
for previously published observations. The ensuing parameters of 
the fit are given in Table[3]and detailed results for each dataset can 
be found in an online-only supplement. The Portalegre transits were 
analysed using the same methods, but only the transit times are used 
below due to the comparatively large scatter of these data. 

The light curves were modelled using the JKTEBOlQ code. 
The primary fitted parameters were the sum and ratio of the frac- 
tional radii of the star and planet, rA+rb and k = and the 
orbital inclination, i. The fractional radii of the components are de- 
fined as TA ~ and = where a is the orbital semimajor 
axis, and Ra and Rh are the true radii of the two objects. Addi- 
tional parameters of the fit included the magnitude level outside 
transit and the midpoint of the transit. 

We generated solutions with each of five limb darkening (LD) 
laws (linear, quadratic, square-root, logarithmic and cubic), and 
with three different treatments of the LD coefficients. The first pos- 
sibility is to fix both coefficients to values predicted using model 
atmospheres; this leads to a dependence on stellar theory as well 
as slightly worse fits due to the larger number of degrees of free- 
dom. The second option is to fit for both coefficients, but this is 
possible only when the data are of extremely high quality. Unless 
otherwise stated, we go for a third alternative: fit for the linear LD 
coefficients and fix the nonlinear one t o theoretically pred icted val- 
ues ('LD-fit/fix' in the nomenclature of Southworthl20 1 C ). The two 



coefficients are highly correlated (e.g. lSouthworth et alj 20073) so 
the theoretical dependence inherent in this approach is negligible. 

Uncertainties in each solution were calcul ated in two ways: 
from 1000 Monte Carlo (MC) simulations jSouthworth et all 
|2004). and with a residual-permutation (RP) algorithm 
jjenkins et al.l l2002h . The larger of the two possible errorbars 
was retained for each fitted parameter. Orbital eccentricity (e) and 
periastron longitude (uu) were incorporated using the constraints 
ecosa; = -0.0099 ± 0.0036 and esino; = -0.0060 ± 0.0069 
dWinn et alj boiObI) . These constraints were treated as obser- 



^ JKTEBOP is written in FORTRAN77 and the source code is available at 
|http : / / www .astro . keele .ac.uk/~jkt/, 



vational data and e cos u! and e sin cj were included as fitted 
parameters. 

3.1 Analysis of each dataset 

The two Cassini transits were modelled together, after scaling 
the observational errors for each transit to give a reduced of 

= 1-0. This step is necessary because the errors returned by 
the aperture photometry routine we use are usually too small. Due 
to the possibility of TTVs in the HAT-P-13 system one has to be 
careful when combining data. In this case we fitted for the orbital 
period (Porb) and the midpoint of the first transit, which is equiva- 
lent to fitting for the two transit midpoints. The RP errorbars were 
selected because they are larger than the MC ones. The LD-fit/fix 
results were adopted. A final value for each photometric parameter 
was obtained by taking the weighted mean of the four values from 
the fits for the non-linear LD laws. Its errorbar was taken to be the 
largest of the four alternative values, with a contribution added in 
quadrature to account for any dependence of the parameter value 
on t he choice of LD law . 

iBakos et alj ( |2009|) presented i-band data obtained with the 

1.2 m telescope and KeplerCam at the F. L. Whipple Observatory 
(FLWO). The 3719 datapoints extend over seven transits within an 
interval of approximately one year, but only two of these transits 
have full phase coverage. We therefore converted the timestamgs 
into o rbital phase (using the ephemeris calculated by Fulton et alj 
l20Tll) . sorted them and combined each set of eight consecutive 
points, to obtain 466 phase-binned points. This will wash out any 
TTVs present ove r that time interval, but inspection of fig. 7 in 
iBakos etal.l ( l2009l) shows that no significant variations exist. A pre- 
liminary fit returned ~ 7.53 so the errorbars were scaled up by 
\/7.53. The LD-fit/fix results are adopted and combined as above. 
The RP errors wer e larger than the M C ones, which is unusual for 
phase-binned data ( ISouthworthl201 ll) but accounted for in our anal- 
ysis. 

ISzabo etai] ( l2010h obtained V- and _R-band coverage of two 
transits; we did not consider the l/-band data of the second transit 
as it suffers from systematic errors. The i?-band observations were 
solved with Porb as a fitted parameter, so the solutions are not sen- 
sitive to the effects of putative TTVs. Given the limited quantity of 
the data we did not attempt LD-fitted solutions. In both cases we 
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found that correlated noise was unimportant and adopted tlie LD- 
fit/fix solutions. The To values from the two datasets (referenced 
to cycle —12) are in poor agreement with each other (see Tables 
A3 and A4) , and r oughly bracket the time of midpoint derived by 
ISzabo et"al] feOlOh . All three timings deviate from the ephemeris 
deri ved below by mu ch more than the derived uncertainties. 

Pal et al.l ^201 1[) presented photometry of three transits from 
two telescopes located at Konkoly Observatory. The first transit was 
obtained with the Schmidt telescope and a CCD camera equipped 
with a Bessell / filter. The second and third came from the 1.0 m 
telescope with VersArray CCD camera and a Cousins R filter. The 
two datasets were modelled separately after scaling up their error- 
bars to enforce Xv = 1-0. For both bands we adopted the LD-fit/fix 
solutions. The MC errorbars are larger than the RP ones, indicating 
that red noise is not important; 

The data presented by iNascimbeni et al. I l l2011bl) were ob- 
tained with the Asiago 1.82 m telescope and AFOSC imager, 
through a Cousins R filter. They comprise 12 585 observations cov- 
ering five closely adjacent transits with cadences ranging from 5.8 s 
to 9.7 s. The data were taken over only 38 days, so can be combined 
without suffering smearing effects due to TTVs with periodicities 
above several months. We therefore phase-binned them by a factor 
of 25 into 504 bins, during which process a 4a clip removed 1 1 of 
the observations. A preliminary fit returned Xv = 1-07 so the sup- 
plied observational errors were left unmodified. The results show a 
strong preference for weaker LD than theoretical expectations, and 
the values for rA and rb depend somewhat on the treatment of LD. 
The LD-fixed solutions can be rejected due to a lower quality of 
fit, and the LD-fitted equivalents return unphysical coefficients, so 
the LD-fit/fix alternatives were adopted. RP errors are smaller than 
MC ones, as we o r dinari ly find for phase-binned data. 

iFulton et al.l ( 1201 ll) obtained photometry of ten transits of 
HAT-P-13, of which only four were fully covered. These include 
two transits taken just over a year apart using the Faulkes Telescope 
North (FTN) and a Z filter, which were solved with a free Porb to 
allow for the possibility of TTVs. A third transit was obtained in the 
i band with the FLWO 1.2 m. Finally, one i-band transit from the 
Sedgwick 0.8 m is accompanied by a partial transit obtained only 
three nights earlier, allowing these data to be modelled with Porb 
fixed to any reasonable value. We found that the errorbars of the 
Sedgwick and FTN data had to be multiplied by \/8.05 and \/5.35, 
respectively, to obtain xS = l-O. For the FTN and Sedgwick data 
we were able to adopt the LD-fit/fix solutions, but for FLWO had 
to stick to LD-fixed as attempts to fit for LD coefficients returned 
unphysical values. For the FTN and FLWO data the RP errors are 
moderately larger than those from the MC algorithm, implying that 
correlated noise is significant in these light curves. 



Table 4. Times of minimum light of HAT-P-13 and their residuals versus 
the ephemeris derived in this work. 



Time of minimum Cycle Residual Reference 

BJD(TDB) - 2400000 no. (JD) 



54581. 


,62443 


±0.00122 


-204, 


.0 


-0.00183 


Bakos et al. (2009) 


54777, 


.01324 


±0.00100 


-137, 


.0 


-0.00099 


Bakos et al. (2009} 


54779, 


.92990 


± 0.00063 


-136 


.0 


-0.00057 


Bakos et al. (2009^ 


54782, 


.84394 


±0.00155 


-135, 


.0 


-0.00277 


Bakos et al. (2009) 


54849, 


.92099 


± 0.00075 


-112, 


.0 


0.00080 


Bakos et al. (2009_) 


54882, 


.00078 


±0.00150 


-101, 


.0 


0.00197 


Bakos et al. (2009) 


54960, 


.74005 


±0.00178 


-74, 


.0 


0.00281 


Bakos et al. ('20091 


55167, 


.79647 


± 0.00280 


-3, 


.0 


0.00631 


Gary(AXA) 


55194, 


.03566 


± 0.00229 


6, 


.0 


-0.00065 


Fulton etal. 120111 


55196, 


.95450 


±0.00127 


7, 


.0 


0.00195 


Fulton etal. ('2011) 


55199, 


.86837 


±0.00123 


8, 


.0 


-0.00041 


Tieman(TRESCA) 


55199, 


.86867 


±0.00131 


8, 


.0 


-0.00011 


Fulton etal. (2011) 


55231, 


.94542 


±0.00091 


19, 


.0 


-0.00199 


Fulton et al. (2011) 


55249, 


.45117 


± 0.00200 


25, 


.0 


0.00634 


Szabo et al. (20101 


55269, 


.86567 


±0.00180 


32, 


.0 


0.00717 


Gary(AXA) 


55272, 


.77577 


±0.00120 


33, 


.0 


0.00103 


Gary(AXA) 


55272, 


.77627 


± 0.00250 


33, 


.0 


0.00153 


Foote(AXA) 


55275, 


.69207 


±0.00180 


34, 


.0 


0.00109 


Gary(AXA) 


55275, 


.69312 


± 0.00266 


34, 


.0 


0.00214 


Fulton etaf (2011) 


55307, 


.77077 


± 0.00370 


45, 


.0 


0.00117 


Gary(AXA) 


55310, 


.69197 


±0.00250 


46, 


.0 


0.00613 


Gary(AXA) 


55511. 


.90854 


±0.00141 


115, 


.0 


0.00226 


Fulton etal. (2011) 


55558, 


.56302 


± 0.00098 


131, 


.0 


-0.00307 


Pal et al. (20JT) 


55561 


.48416 


± 0.00400 


132, 


.0 


0.00183 


Pal et al. (2011) 


55564, 


.39876 


±0.00180 


133, 


.0 


0.00019 


Nascimbeni et al. (201 lb) 


55584, 


.81455 


±0.00153 


140, 


.0 


0.00231 


Dvorak(TRESCA) 


55590, 


.64523 


±0.00179 


142, 


.0 


0.00052 


Pal etaf (2011) 


55593, 


.55879 


±0.00185 


143, 


.0 


-0.00216 


This work (Portalegre) 


55593 


.56147 


+ 001 1 5 


143 


.0 


00052 


NaRPirnhpni pf al f^901 1 


55596, 


.47291 


±0.00140 


144, 


.0 


-0.00428 


Naves(TRESCA) 


55596, 


.47327 


± 0.00202 


144. 


,0 


-0.00392 


This work (Portalegre) 


55596. 


,47662 


±0.00305 


144. 


.0 


-0.00057 


Nascimbeni et al. (201 lb) 
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This work (Cassini) 
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Fulton etal. (20111 


55622. 
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±0.00166 


153. 


,0 


0.00018 


Fulton etal. (20111 


55669. 


,38140 


±0.00126 


169. 


,0 


-0.00175 


This work (Cassini) 



3.2 Combined results 

The photometric parameters resulting from the ten light curves are 
collected in Table[3l We calculated the weighted means to obtain 
final values of these quantities. The Xv values of the averaging pro- 
cess are all below 0.5, with the exception of k. We deduce that the 
ten light curves are in sufficient mutual agreement. The final values 
are dominated by the results from the Asiago light curve, which is 
easily the best of the available datasets. 

For k the averaging process yields Xu = 2.0, so the agree- 
ment between the light curve solutions is not so good. The Cassini 
data are the primary contributor to this situation, as they point to a 
larger k than the rest of the datasets. Individual solutions of the two 
Cassini light curves yield similar values of k. Moderate disagree- 



ments in k have been frequently seen in the Homogeneous Studies 
papers and can be attributed to starspot activity and/or systematic 
errors in the light curve, k is also determined to a high precision 
(by comparison to rA, r^, and i) so systematic differences are com- 
paratively obvious. 

Our final photometric parameters (Table|3]l are so mewhat dif- 
ferent to values based o n the discovery photometry dSakos et al.l 
2009, ; IWinnetalj2010bl) . In particular we find a lower i and larger 
rA- The latter quantity is observationally strongly tied to rb so 
our results point towards a larger star and planet than previously 
proposed. Our photometric parameters agree re asonably well with 
those proposed recently bv iFulton et a l|), but have smaller 
uncertainties. 
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Figure 3. Plot of the residuals of the timings of mid-transit o f HAT-P- 1 3 vers us a linear ephemeris. The timings in b lack are fr om this work , in gre y are from 
iBakos et alj |2Q09 ). blue from lSzabo et al.l I201Q) . hlac from lPa et al] j201ll) . green from Nascimbeni et al.l i2011bl) . r ed fromlFulton et alj j201lh . and open 
circles for the amateur timings. The solid line shows the ephemeris fro m the current work and the dotted line that from lBakos et alj j2009l) . The dashed curve 
is an approximate representation of the possible periodicity proposed by lNascimbeni et al] j2011lj) . 



4 ON THE TRANSIT TIMINGS OF IIAT-P-13 B 

The third body in the HAT-P-13 system, known to have an eccen- 
tric orbit with = 0.6616 ± 0.0052 and Po^.c = 446.22 ± 
0.27 d from radial velocity measurements, is expected to cause 
TTVs within the in n er HAT-P-13 A,b system ( iBakos et alj|2009l : 
IPavne & Fordl20l"lh . lBakos et alj found "suggestive" but not "sig- 
nificant" evidence for TTVs in t heir data obtained i n their 2007- 
8 and 2008-9 observing seasons. ISzabo et"ail ( 1201 Oh obtained an 
improved orbital ephemeris and placed an upper limit of 0.001 d 
on size of the phenomenon, with the inclusion of their two transit 

times from the 2009-10 season. 

However, the times of three transits acquired by |Pal et al.l 

J) during the 2010-11 season were about 22 min utes early 
with respect to the previous orbital ephemerides, leading Pal et al.| 
to claim a significant detec tion of TTVs. The timin gs found by 
|Pal et alj were supported by Nasci mbeni et al.l ( 1201 l b) on the ba- 
sis of five new tra nsits obtained during the same observing season. 
iNascimbeni et al.l demonstrated that a sinusoidal TTV function pro- 
vide d a good fit to all ex isting timing measurements. 

iFulton et al.1 ( 1201 ih have subsequently presented ten transit 
timings which cast doubt on the possibility of TTVs: five from 
the 20 10-11 season which ag r ee wel l with those from |Pal et al.l 

If) and INascimbeni et al T bOllbh . and five fror n the previ- 
ous se a son which co nflict with the timings found by ISzabo et al.l 
( l2010h .l Fulton et alj reanalysed all published follow-up observa- 
tions of HAT-P-13 and concluded that they were consistent with a 
linear ephemeris w ith the exception of the first transit dataset from 
ISzabo et aljfcoid) . 

In order to firmly establish the character of the situation, we 
have collected all available transit midpoint times for the HAT-P- 
13 A,b system. We have used the timings as quoted by the orig- 
inal sources , rathe r than adopting those from the reanalysis by 
iFulton et al. U20Tlh . We included ten timings obtained by amateur 
astronomers and placed on the AX^fl and TRESC^] websites. We 



^ Amateur Exoplanet Archive,lhttp : / /brucegary . net/AXA/ x . htm 
^ The TRansiting ExoplanetS and CAndidates (TRESCA) website can be 
found at, http : / /var2 . astro. cz/EN/tresca/ index, php 



rejected amateur timings which are based on data that are either 
very scattered or do not cover a full transit. 

Transit timings were obtained from our own observations by 
fitting JKTEBOP models to each transit following the LD-fit/fix pre- 
scription. The flux normalisation was allowed to vary linearly with 
time. Errors were estimated from RP and from 1000 MC simula- 
tions, and were multiplied by two in order to guard against any 
undetected systematic noise in the data. At the request of the ref- 
eree we also asse ssed correlated noise using the '/?' approach (e.g. 
IWinn et al. 1(2007). We evaluated values for individual transits and 
for groups of between two and ten datapoints, finding a maximum 
/3 of 1.26. The corresponding increases in the uncertainties in the 
Tq values are smaller than the factor of two we used above. 

All timings were placed on the BJD(TDB) time system. We 
fitted a straight line to obtain a new orbital ephemeris, finding 
xi = 2.00 with one obvious outlier. After the rejection of th e 
offending point, which is the first timing from lSzabo et alj ( l2010h . 
we obtained: 

To = BJD(TDB) 2 455176.53878(27) + 2.9162383(22) x E 

with Xu = 1.54. The bracketed quantities represent the uncertain- 
ties in the ephemeris, and have been increased to account for the 
excess Xv- The full list of transit timings and their references is 
given in Table|4] It should be noted that many of the timings are 
measured from data covering only part of a tran sit, which is known 
to reduce their reliability (e. g.lGibson et al. 2009). 

The two timings from lSzabo et al.. C2010;) both deviate from 
the orbital ephemeris above, being later by 8.6a and 3.2(7. Our own 
analyses of these data return timings which are similarly distant 
from expectat ions. A detai l ed rea nalysis of the corresponding data 
performed b V IFulton et"^ bOllh resulted in a timing for the sec- 
ond of these transits which conflicts less with a linear ephemeris 
(1.6(t). The discrepancy of the first transit remains unexplained. A 
few of the amateur timings are late by a similar amount to this one, 
but with much larger errorbars. 

We conclude that the available data do not provide a clear in- 
dication of the existence of TTVs, primarily on the basis that we 
cannot conceive of a reasonable TTV function which is a signifi- 
cant improvement over a linear ephemeris. The transits which oc- 
cur later than predicted by our ephemeris are not grouped together. 
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Table 5. Final physical properties of the HAT-P- 1 3 system, compared with results from the literature. Where two eiTorbars are given, the first refers to the 
statistical uncertainties and the second to the systematic errors. 



This work (final) Bakos et al. (2009) Winn et al. (2010b) 



Ma(Mq) 


1.320 ± 0.048 ±0.039 


1 r,iQ + 0.050 

— 0.099 


1 99+0.05 
^■Z'^ — 0.10 


i?A(R0) 


1.756 ± 0.043 ±0.017 


1.559 ±0.082 


1.559 ±0.080 


log OA (cgs) 


4.070 ± 0.020 ± 0.004 


4.13 ± 0.04 




PA (Pq) 


0.244 ±0.013 






Mb ( Mjup) 


0.906 ± 0.024 ±0.018 


n ac;q+0.029 


0.851 ±0.038 


Rb (Rjup) 


1.487 ± 0.038 ±0.015 


1.281 ± 0.079 


1.272 ±0.065 


9b (ms~i) 


10.15 ±0.43 


12.9 ± 1.5 




Pb (P,Tup) 


0.257 ± 0.017 ±0.003 


„ q7c:+0.078 




Te'q (K) 


1725 ± 31 


1653 ± 45 




e 


0.0404 ± 0.0023 ± 0.0004 


0.046 ± 0.003 




a(AU) 


0.04383 ± 0.00053 ± 0.00043 


0427+°-°°°'^ 




Age (Gyr) 


q +0.3 

•J-'J -2.9 -0.7 


5 n+2-5 





but are interleaved with ones which happen at the expected times. 
An explanation involving TTVs therefore would require a highly 
contrived functional form. 

S o where did th e previous suggestions of TTVs come 
from? lPaletal.1 (1201 ih used an earlier ep hemeris, tuned on the 
iBakos et alj ( l2009h and lSzabo etZH201(]i) observations, to show 
that their transit timings were earlier than expected. The dotted 
line in Fig.[3] represents this ephemeris and shows that it fails to 
match the more recent transit timings. Nascimbeni et al. ( l2011bl) 
suggested that a sinusoidal TTV of amplitude 0.005 d and pe- 
riod 1150d was in good correspondance with the observations, as 
demonstrated by their fig. 2. We have endeavoured to place this pe- 
riodic variation, whose parameters were not fully specified, onto 
Fig.|3] with a little manual fine- tuning. The 2009-10 transit timings 



obtained bv lFulton et al.l (1201 1 ) clearly dismiss the sinusoidal TTV 
proposed by Nascimbeni et al. i2011bh . leaving a linear ephemeris 



as the only reasonable option. 



5 PHYSICAL PROPERTIES OF THE HAT-P- 13 SYSTEM 

Several sets of information exist from which the properties of 
the HAT-P-13A,b system can be derived. Analysis of the avail- 
able light curves has given values for P prb, i, ta and rb- The 
high-precis i on radi al velocities procured bv lBakos et al. |(|2009') and 
IWinnetal.ll2010bl) supply measurements of the velocity amplitude 
of the star {Ka = 106.04 ± 0.73 ms~^) and the orbit al eccentric- 
ity (e = 0.0133 ± 0.0041). Analysis of the spectra bv lBakos et"al] 
( l2009l) furthermore lead to estimates of the stellar effective temper- 
ature (Tch = 5653±90K)andmetallicity([Sl] = +0.41±0.08). 
Finally, constraints on the properties of the star can be obtained by 
interpolation within tabulated predictions from stellar evolutionary 

models. 

Our solution process ( ISouthwortlj|2009l) consists of finding 
the best agreement between the observed and model-predicted 
TeffS, and the measured ta and calculated -5^. This is done us- 
ing the velocity amplitude of the planet, K^, as a solution control 
parameter and calculating the full system properties using stan- 
dard formulae (e.g. .Hild itch 2001.1 . The system properties com- 
prise the mass, radius, surface gravity and mean density for the 
star (Ma, Ra, log pa and pa) and planet (Afb, ^^b, Pb and pb), 
the orbit al semima j or axi s (a), the planetary equilibrium tempera- 
ture and ' Safronovl ( [1973) number {T^^, O) and an estimate of the 
evolutionary age of the star. 



The statistical errors on the resulting values are calculated us- 
ing a perturbation analysis dSouthworth et alj|2005h which yields 
a full error budget for each output quantity. The use of theoretical 
models incurs a dependence on stellar theory which is assessed by 
running separat e solutions with e ach of five different sets of model 
tabulations (see ISouthworthll20rQ for details). Finally, an alterna- 
tive empirical es timate of the phys ical properties is obtained us- 
ing a calibration dSouthworth" 200 ) base d on eclipsing binary star 
systems, inspired by Enoch et al] 1 201(]|) and|^rres et al. (2010). 
The control parameter A'b represents the constraints obtained from 
stellar theory. Note that the values ffb, Pa and T^^ are not re- 
liant on constraints from stellar theory ([Southworth et al.ll2007bl ; 
ISeager & Mallen-Ornelasll2003l : ISouthwortriboTll) ! 

The sets of physical properties arising from each of the five 
stellar model tabulations and from the empirical calibration are 
given in Table All. For the final results in Table|5]we adopted the 
unweighted mean of each parameter from the solutions for the five 
stellar models. The statistical error is the largest of the individal er- 
rorbars, and the systematic error is the standard deviation of the 
five values. Compared to previous studies, we find a larger and 
more evolved star and a correspondingly slightly more massive but 
significantly bigger and hotter planet. The extensive photometric 
dataset considered in this work leads to more precise measurements 
of the physical properties, but the uncertainties in ta and rb con- 
tinue to dominate the error budget. The uncertainty in A/a stems 
mainly from that in , suggesting that a new spectral synthesis 
study would also be useful in improving our understanding of the 
HAT-P- 13 system. 



6 SUMMARY 

The HAT-P- 13 system is unusual in that a transiting hot Jupiter and 
its host star are accompanied by a clearly-detected third compo- 
nent on a wider orbit. Such a configuration should result in HAT-P- 
13 c inducing TTVs within the HAT-P-13 A,b system, which may 
be detectable within a comparatively short time period. This pos- 
sibility has generated substantial interest, resulting in a large body 
of photometric observations covering many transits of the star by 
the inner planet. We have assembled the available transit timing 
measurements and shown that they are most easily explained by a 
linear ephemeris, albeit with a small number of values which occur 
later than expected. The discrepant measurements are not clumped 
together, so could only be explained via highly complex functional 
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forms. Previous claims of TTVs can be attributed to small-number 
statistics, although continued photometric monitoring has a good 
chance of turning up something interesting in the future. 

We have presented new observations of four transits, obtained 
using telescope defocussing techniques. Including previously pub- 
lished data, we have ten good sets of transit light curves. These 
were each analysed within the cont ext of our H omo/^eneoiis Stud- 
ies project ( Southworth 2008. 2009 „201(]|.l201 Ih . and a good agree- 
ment between the results was found. We combined them with the 
measured spectroscopic properties of the host star and several sets 
of theoretical stellar model predictions to find the physical proper- 
ties of the system. HAT-P-13 is now well-characterised, although 
additional photometric and spectroscopic measurements would al- 
low further improvement. We have included it in the TEPCat cata- 
logu^fl of the physical properties of transiting planetary systems. 

We find a significantly different set of physical properties 
compared to previous studies, which had access to only two of the 
ten photometric datasets used here. The star is more massive, larger 
and more evolved. The planet, whose properties are measured rel- 
ative to its host star, is similarly heavier and bigger. Its lower den- 
sity and higher equ ilibrium temperature place it firmly in the 'pM' 
class advocated bv lPortnev et al.l l l2008h . Its radius is too large to 
m atch the values pred icted by the models of iFortnev et alj j2007l) 
or lBaraffeet"ai](l2008h. 

iLaughlin et al. f Hoilh found that the radius anomaly (the mea- 
sured radius of a TEP versus that predicted by theoretical mod- 
els) is correlated with equilibrium temperature, and possibly in- 
versely correlated with host star [^] . The large radius anomaly 
and high equilibrium temperature of HAT-P-13 b corroborate the 
former observation, but the highly metal-rich nature of the parent 
star ( [^] = 0.41 ± 0.08) is contrary to the latter suggestion. 
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